Intrinsic instability of electronic interfaces with strong Rashba coupling 
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We consider a model for the two-dimensional electron gas formed at the interface of oxide het- 
erostructures, which includes a Rashba spin-orbit coupling proportional to the electric field per- 
pendicular to the interface. Based on the standard mechanism of polarity catastrophe, we assume 
that the electric field has a contribution proportional to the electron density. Under these simple 
and general assumptions, we show that a phase separation instability (signaled by a negative com- 
pressibility) occurs for realistic values of the spin-orbit coupling and of the electronic band-structure 
parameters. This provides an intrinsic mechanism for the inhomogeneous phases observed at the 
LaAlOa/SrTiOa or LaTiOa/SrTiOa interfaces. 

PACS numbers: 71.70.Ej,73.20.-r,73.43.Nq,74.81.-g 
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The observation of a two-dimensional (2D) metal- 
lic state at the heterointerface of two insulating oxides 
[LaAlOa/SrTiOg (LAO/STO)] [l| remarkably discovered 
a new class of high-mobility electron gases (EGs) impor- 
tant both for fundamental and applicative reasons. The 
occurrence of superconductivity in this 2DEG 0, Q and 
in LaTiOa/SrTiOa (LTO/STO) H with the possi- 
bility to tune the charge density by gating, has further 
attracted great attention. On the other hand, there is 
increasing evidence that electron inhomogeneity plays a 
relevant role in these systems. Not only the large width 
of the superconducting transition in transport experi- 
ments is a clear indication of charge inhomogeneity 
but also magnetometry and tunnelling experiments 
JjJ find submicrometric inhomogeneities. While impu- 
rities, defects, and other extrinsic mechanisms (see, e.g. 
Ref. [l3|) surely induce inhomogeneities, the recent dis- 
covery of negative compressibility in a low filling regime 
[l3| is a stringent demonstration that an intrinsic mech- 
anism (like, e.g., an effective attraction) is at work to 
render these 2DEGs inhomogeneous by charge segrega- 
tion and phase separation even in a perfectly clean and 
homogeneous system. Moreover, even if such mechanism 
were not strong enough to drive the system unstable, it 
would increase the charge susceptibility, emphasizing the 
effects of the extrinsic mechanisms (impurities, defects 
and so on). 

In this Letter, we point out that a generic source of 
phase separation is provided by the Rashba spin-orbit 
coupling (RSOC), whenever the electric field determin- 
ing this coupling also controls the electron density. This 
is precisely the case of LAO/STO and LTO/STO (generi- 
cally, LXO/STO) interfaces, where the two following con- 
ditions remarkably meet: i) a strong electric field occurs 
perpendicular to the interface because of the polarity 
catastrophe 14-1^ and brings the electrons at the ox- 
ide interface to produce the 2DEG. An additional elec- 
tric field, although quantitatively less important, is intro- 
duced by the gating potential and tunes the density of 



the 2DEG; ii) the parity symmetry is broken at the inter- 
face naturally entailing a RSOC, which experiments have 
found to be substantial [l7, 18|. We show that these two 
concomitant conditions are enough to drive the 2DEG 
unstable towards phase separation, thereby providing a 
general and intrinsic mechanism for the inhomogeneity 
of these oxide interfaces. 

The mechanism for RSOC-induced phase separation is 
rather simple. In a metallic system with a rigid band 
structure the chemical potential /i increases upon in- 
creasing the electron density n and the compressibility 
K = dn/dfi is positive. On the other hand, if the band 
structure is modified by the charge density (like, e.g., 
in strongly correlated systems, where the quasiparticle 
bandwidth increases moving away from the half-filled 
Mott-Hubbard insulator) the possibility may occur that 
/i decreases with increasing n and k < 0, with the nega- 
tive compressibility region signaling a phase separation. 
The mechanism at work is so generic that it could be rel- 
evant for many other systems in the presence of RSOC 
like, e.g., surface states of topological insulators, InAs 
or GaAs semiconductor heterostructures, or heavy metal 
surface alloys. 

Our analysis finds that a band structure with large 
(even anisotropic) masses together with a rather isotropic 
Rashba coupling are favorable conditions for this insta- 
bility to occur. We also stress that the charge inhomo- 
geneity within the interface is due to an inhomogeneous 
electronic reconstruction and therefore it is balanced by 
a correspondingly inhomogeneous redistribution of the 
countercharges on the top layer. Therefore the inhomo- 
geneity is not prevented by the standard Coulomb mecha- 
nisms leading to frustrated phase separation in other sys- 
tems 3" 21 1 . This naturally leads to large submicromet- 
ric inhomogeneities, like those detected in LXO/STO. 

Fig. [TJa) schematically shows the 2DEG between the 
two oxide layers, the gating layers and the resulting elec- 
tric fields. The electron reconstruction to avoid the polar- 
ity catastrophe is also depicted with the electron charge 
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FIG. 1. (Color online) (a) Schematic view of the LXO/STO 
interface in the presence of gate potential Vg. no is the density 
of the electrons transferred to the interface by the polarity 
catastrophe and giving rise the Epc electric field. Sn is the 
electron density tuned by the gating field Eg. (b) Schematic 
band structure for the t2g bands along with = 30 and 
rj — 1 and (c) rj — l/i^. We report the d^y (solid black), d^z 
(dashed blue online) , dy^ (dot-dashed, red online) bands. 



density no transferred onto the interface, leaving oppo- 
sitely charged planes in the LXO layers above it. As a 
consequence, a field Epc = eno/(eoei) arises perpendic- 
ular to the interface, where ei ~ 20 is the LXO dielectric 
constant, eg is the vacuum dielectric permittivity and e 
is the electron charge. Typical values of Epc may reach 
lO^V/m. The total electric field perpendicular to the 
interface is obtained by adding to Epc the gate elec- 
tric field Eg = Vg/d, where d is the width of the STO 
substrate, which tunes the charge density with a term 
eSn = eQ€2Eg. For a STO substrate 0.5 x 10^'^ m thick. 
Eg ~ 10^ V/ni (<C Epc)- It is quite important to notice 
that the large value of the low-temperature STO dielec- 
tric constant £2 ^ 10^ allows for relatively large varia- 
tions of the electron filling 5n, while the gate electric field 
always stays much smaller than Epc ■ As a consequence, 
the 2DEG can be greatly depleted when Sn ~ —uq, while 
the total electric field is still large E = Eg + Epc ~ Epc- 

For a description of the STO band structure, where the 
electrons of the LXO/STO interfaces mainly reside, we 
consider three parabolic bands schematizing the bottom 
of the t2g-like bands of Ti. The lowest (c?xj/) band has 
light isotropic mass ~ 0.7mo (mg is the bare electron 



mass) while at an energy A « 50 meV above there are 
two {dxz and dyz) bands with the same light mass in one 
direction and a heavy mass as large as 20mo in the other 
[isl . [i^ - E?! ] . For simplicity, we neglect all mixings among 
the three bands 

The lowest band has therefore isotropic dispersion 
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which is the prototypical model for 2DEGs with Rashba 
coupling a (see, e.g., Ref. fso']) and gives rise to a 2D 
isotropic band structure composed of two branches [solid 
black curves in Fig. [T] (b),(c)] split by — 2a\k\, 
with a minimum —Eq = —ma^/{2h^) occurring on a cir- 
cle of radius am. The corresponding density of states 
(DOS ) is N{e) = Nq = m/inh^) for e > and N{e) = 
TVq/a/I + e/so for — £0 < e < 0. The other two bands 
give rise to (strongly) anisotropic dispersions 
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with V = my/mx ^ 30 and 77 = ay /ax for the dxz band 
and = mx/my ^ 30 and 77 = ax/ay for the dyz band. 
The DOS arising from Eq. ([2]) depends on (e — A)/£o 
only and is given in terms of complete elliptic integrals 
of the first and third kind (see Supplemental Material). 
The precise form of ax and ay, and therefore of rj, should 
be determined by first-principle calculations, which are 
beyond the scope of this Letter. However, we argue that t] 
should be intermediate between two extreme cases: rj = 1 
(isotropic Rashba coupling) and rj = 1/v [i.e. ax,y = 
^/inx,y to reconstruct the relativistic form of the spin- 
orbit coupling: (v x cr) • E] . The band bottom occurs at 
A — vEo in the former case and at A — eo in the latter 
(see Fig.mb) and (c)). 

In any case, a and ax,y depend on the perpendicu- 
lar electric field E but, instead of the customary sim- 
ple proportionality, we here take a non-linear expression 
a{E) = aE/ (1 -I- /3E)^ which is also valid at large fields 
and stems from a standard derivation of the Rashba 
couphng [sO, im (see also, e.g., Eq. (2) in Ref. [sl ] 
and Supplemental Material of this Letter). Thus, £0 = 
'^E'^/ {1 + /3Ef [with 7 = mce'^ / {2h'^)] is an increasing 
function of E up to moderate-large values [~ 1/(2/?)]. 
The key point here is that the electric field is directly 
related to the number of electrons in the plane E{n) = 
E{0) + nE' and therefore the larger is n, the deeper is the 
energy minimum £0 ~ [E{0)+nE']'^ . This may render en- 
ergetically convenient for the system to attract electrons 
to have a deeper energy minimum, where more electrons 
can be accommodated at lower energy. This downward 
shift of the band bottom may overcome the increase of 
the Fermi level due to the increased n and an overall 
decrease of may occur, leading to a negative compress- 
ibility. 



Filling the bands in Eqs. (Ill2p one easily gets the gen- 
eral expression of /it as a function of n = ng + 5n at a 
fixed gate potential Vg. At lower densities, when only 
the isotropic band is involved, one obtains 

— -2eoKl^3) M>0. 

Experiments are usually performed at a fixed Vg (and 
therefore at fixed Sn) and a stable system uniformly re- 
constructs the surface forming a uniform uq electron den- 
sity at the interface to avoid the polarity catastrophe. We 
now ask whether the system could be unstable and, at 
fixed Vg, could display a tendency to inhomogeneous re- 
constructions and a negative compressibility upon vary- 
ing n (or, equivalently, uq since dn is kept fixed). 

The inverse compressibility k^^ = dfi/dn at fixed Vg 
for the isotropic band reads 
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and can be negative in both the low-density (/i < 0) 
and in the high-density {fi > 0) regime. The condition 
K < naturally implies an electronic phase separation, 
i.e. a separation of the system into regions of differ- 
ent electronic densities to be determined by the standard 
Maxwell construction on the fi vs. n curve. 

An analysis of Eq. (4) for /i < shows that if the 
electric field (and therefore a) stays finite when n = no + 
Sn (i.e. E{0) ^ 0), then k"! -deo/dn < 
and the system is always driven unstable at sufficiently 
low filling. Although, depending on the value of a, this 
latter instability could occur at very low densities [the 
Sn — —0.02 curve (red online) in Fig. HKb) displays a 
barely visible negative initial slope around rip = 0.02 (i.e. 
n « 0)] , this observation might turn out to be relevant in 
physical systems like, e.g., MOSFETs or semiconducting 
heterostructures (see, e.g., Ref. [s^]) where a is small 
but n also is very small. In this regard, we considered 
the possibility that RSOC mechanism could be at the 
origin of the negative compressibility recently observed 
in LAO/STO interfaces ^ at low filling. We find that a 
negative compressibility in the observed range n ^ 10^^ 
cm~^ would require values of a ^ 5 x 10^^^ eVm. 

On the other hand, in the case > of Eq. (4) we 
can take the limit of not too large fields f3E <^ 1 [and 
E « nE'] and obtain k^^ « n~^(n/iVo — 4eo) which 
gives the simple condition n/No < Aeq for the negative 
compressibility to occur. For instance, at densities n ^ 
10^^ cm~^, K < would require values as large as a ~ 
6 X 10~^^ eVm. It is also important to notice that a large 
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FIG. 2. (Color online) Chemical potential as a function of 
reconstructed electron density no [0.01 electrons/cell corre- 
spond to 6, 25x 10^^ cm-2] for the three t2g bands of STO. The 
lowest band is isotropic {u = = 1) with mass m ~ 0.7mo. 
The two bands lying A = 50 meV above are anisotropic with 
v — 30 and rj = 1 {a) or rj = l/v (b). The legenda report 
the values of the additional electrons induced by gating. In 
(a) a Maxwell construction is also reported for the Vg = 
case. The dotted curves in (a) and (b) (right axes) are a as 
a function of n = no (Vg = 0). 



DOS Nq decreases the positive contribution to k ^ and 
favors the occurrence of the instability. 

Recent magneto-conductivity experiments in 
LAO/STO interfaces find substantial values for 
a ~ 10-12 10-"eVm [U IH. Similar values 
have been obtained for LTO/STO interfaces [s^. These 
values are about 5-8 times smaller than those estimated 
above to drive the system unstable. In the specific case 
of very small densities [l^j , this quantitative discrepancy 
might well be due to the lack of many-body effects and 
disorder in our schematic model 3J| ) . The question then 

2-=-4x 10^3 cm-2, 



arises whether at higher densities n 
with the anisotropic bands partially filled, LXO/STO 
interfaces displays a negative compressibility at lower 
values of a. This question is addressed by Fig. [51 which is 
the main outcome of our analysis and demonstrates the 
possibility of an inhomogenous electronic reconstruction 
as soon as the anisotropic bands begin to fill. While 
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the rj = 1/v case [panel (b)] only becomes unstable 
for large values of a ~ 1 5 x 10~^^eVm, the rj = 1 
case [panel (a)] displays electronic phase separation for 
quite realistic values of a ~ 0.1 -;- 1 x 10~^^ eVm, well 
within the range of those found in Refs. [iTI and [iSj in 
LAO/STO. 

To provide a clear understanding of this result, we con- 
sider the following expressions 

n = 2Noeo [1 + 2f{i^, rj)] + fiNo + 2(/i - A)NoV^, (5) 
1 ( n 



A* 



1 + [ iVo 



2eo[l + 2/(i/,ry)] + 2A0; , (6) 



valid for ^ > A and n > 2iVoeo[l + 2/(1^,77)] + AA^o, 
respectively. Here, f(v,ri) is a (rapidly) increasing func- 
tion of the mass anisotropy v (see Supplemental Mate- 
rial) and gives a moderate [/(30, 1/30) ~ 3] to strong 
[/(30,1) ^ 90] contribution to the negative term in 
Eq. (jG]). This shows that the RSOC-mediated instability 
takes advantage from the mass anisotropy and a negative 
compressibility is now more likely to occur. 

For /J, < A the linear growth of in Fig. [5] shows that 
the DOS of the isotropic d^y band is too small to pro- 
mote an instability for the values of a{n) given by the 
dot-dashed curves. On the other hand, as soon as the 
density uq + dn is large enough, the chemical potential 
eventually enters the anisotropic bands, the DOS rapidly 
increases (see, e.g., Ref. [29]) and ^ decreases with no 
signaling a negative compressibility region. For instance, 
the Vg — curve in Fig. [2{a) has negative slope in the 
range hq ~ 0.02 - 0.04 electrons/cell (~ 3 - 5 x 10^^ 
cm~^), which are typical values of as-grown systems. 
Even for lower values of uq, the instability eventually 
occurs upon increasing the density with Vg (see, e.g., the 
= 112 V (blue online) dashed curve). This leads to 
a thermodynamic phase separation into regions of differ- 
ent densities (e.g., ni — 0.022 and n2 = 0.058 in the 
= curve of Fig. Efa)) determined by a standard 
Maxwell construction [36]. For any given filling n, with 
ni < n < n2, the Maxwell construction also determines 
the relative weight of the two phases at ni and n2. 

A few remarks are now in order. Firstly, we point out 
that the band structure of the 2DEG at the LXO/STO 
interfaces is substantially modified by the RSOC arising 
from the strong electric field due to the polarity catas- 
trophe. Therefore, it is not surprising that angle-resolved 
photoemission experiments 24-2^ on STO/ vacuum in- 
terfaces do not show the strong band splitting that one 
would expect from the experiments of Refs. [itI [l8|. 

Secondly, we stress that additional mechanisms are 
also present in the real systems, which may cooperate 
with the intrinsic RSOC mechanism to produce inhomo- 
geneous charge distributions. Indeed, if the electronic 
system has a large (although still positive) compress- 
ibility, impurities and defects more easily induce large 
inhomogeneities in such a "soft" and largely fluctuat- 



ing electron gas. In other words, they act as "external 
fields" on the density and may enhance the effects of the 
RSOC to locally induce phase separation. Previous work 
in strongly correlated systems 21], |37h41|. also shows 



that electron-electron correlations and electron-phonon 
coupling favor phase separation. In this regard, the re- 
cent analysis of RSOC within Hartree-Fock approxima- 
tion [i^l , can be a starting point for investigating the role 
of interactions on the RSOC-mediated instability. 

Finally, we mention that our schematization of the 
LXO/STO interface in terms of a 2DEG should natu- 
rally be relaxed to account for the finite width of the 
interfacial electron gas 23142 7|. However, the introduc- 
tion of subbands due to this finite width will surely in- 
crease the DOS and consequently decrease the /i vs. n 
slope when the chemical potential progressively crosses 
the subbands. This will strongly favor the instability 
further reducing the RSOC required for it. 

In conclusion, within our rather essential description 
of the LXO/STO interface, we find that electronic phase 
separation may likely occur at oxide interfaces, where 
substantial density- dependent electric fields can arise. We 
stress that a negative compressibility (cf. Ref. 113]) can- 
not be explained by impurities, defects and so on, but 
is a distinct signature of an intrinsic mechanism leading 
to an effective charge segregating attraction like the one 
proposed here. We also show that (anisotropic) large 
masses (and subband splitting), isotropic RSOC, and 
low fillings favor the instability. This provides a natu- 
ral framework accounting for the inhomogeneous phases 
observed in LXO/STO interfaces. It would be interesting 
to investigate the possible occurrence of electronic phase 
separation in cases where the RSOC is sizable and/or the 
density is very low, like in quantum wells, in boundaries 
of heavy metal alloys like Bi^^Pbi-a;, in reconstructed sur- 
faces of Ag(lll) 43, 44J, in MOSFET and semiconduct- 
ing heterostructures at low densities, and in surface states 
of topological insulators [s^, ll^ . 
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